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Abstract 

Equations built on fractional derivatives prove to be a powerful tool in the de- 
scription of complex systems when the effects of singularity, fractal supports, and 
long-range dependence play a role. In this paper, we advocate an application of the 
fractional derivative formalism to a fairly general class of critical phenomena when 
the organization of the system near the phase transition point is influenced by a 
competing nonlocal ordering. Fractional modifications of the free energy functional 
at criticality and of the widely known Ginzburg-Landau equation central to the clas- 
sical Landau theory of second-type phase transitions are discussed in some detail. 
An implication of the fractional Ginzburg-Landau equation is a renormalization of 
the transition temperature owing to the nonlocality present. 
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The formulation of fractional kinetics [1] has led to a considerable progress in 
our vision of complex systems at the microscopic level. The use of fractional 
derivative operators — unconventional analytical tools extending the familiar 
high-school calculus to a broader class of mathematical objects — has been 
recognized in the modern physics [2] owing to its elegance and the proximity 
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to the standard, well established algorithms. The fractional kinetic equations 
dealing with generalized derivatives in space and time incorporate in a nat- 
ural, unified way the key features of non-Gaussianity and long-range depen- 
dence that often break down the restrictive assumptions of locality and lack 
of correlations underlying the conventional statistical mechanical paradigm. 
Prom a probabilistic standpoint, fractional kinetics extends Gaussian stochas- 
tic processes (i.e., Brownian random walks) by taking into account long-range 
correlated events in the tail of the probability density function. Such events 
dominate, for instance. Levy-type processes [3] and fractal time random walks 
(FTRW's) [4], the simplest model realizations accounting for anomalous trans- 
port phenomena in turbulent media [5] . Manifestations of fractional kinetics 
have been found in, e.g., point vortex flows [6], low confinement mode plasmas 
[7], non-Gaussianity of fiuctuations measured in the edge and scrape-off layer 
region of fusion devices [8] , etc. A comprehension of the essential role played 
by FTRW and Levy statistics in the microscopic description of turbulence and 
chaos stipulated fractional generalizations of the diffusion and Fokker-Planck- 
Kolmogorov equations, discussed in a series of publications [9,10,11,12,13]. 
Beside the theory of turbulent diffusion, applications of fractional kinetics con- 
cern the fractional Kramers problem [14], relaxation in polymer systems and 
rebinding phenomena in proteins [15] , scale-invariance and universality near a 
phase transition point [16], cosmic rays acceleration [17], the dynamics of frac- 
ton excitations [18], including modulational instability [19] and self-focusing 
of waves on fractals [20], and many other realizations [2]. The current state of 
the art is summarized in review articles [19,21,22]. 

In this Letter, we advocate an application of the fractional derivative formal- 
ism to the thermodynamics of second-type phase transitions in the presence of 
a coexisting nonlocal ordering which may influence the properties of the basic 
"symmetric" phase below the transition point. As an example, we mention 
highly correlated electron liquid states, such as electron liquids with frac- 
tionally charged excitations [23], as well as high-temperature superconducting 
fluid phases in copper-oxide compounds and their derivatives [24] where the 
nonlocal ordering can be associated with the so-called "stripy" order [25]. A 
growing evidence of stripes and of their role in the superconducting transition 
at high temperatures is a hot topic in condensed matter research [26] . 

In what follows, we imply a phase diagram which accommodates two thermo- 
dynamically distinct phases, the "symmetric" (superconducting) phase below 
the transition temperature T^, and "asymmetric" (normal) phase above the T^. 
As usual, the deviation between the symmetric and normal phases for T — > Tc 
is characterized by the order parameter, ip — il){x), which is assumed to be an 
analytical function of a 1-dimensional coordinate variable, x. 

Our study has two interconnected goals. First, we suggest that an interaction 
between the ip order and a coexisting nonlocal ordering can be characterized 
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by a fractional generalization of the free energy expansion near the transition 
point. Second, we demonstrate that the order parameter ip obeys an integro- 
differential equation which can be fairly considered as a fractional extension 
of the widely known Ginzburg-Landau equation, central to the classical Lan- 
dau theory of second-type phase transitions. This unconventional, fractional 
equation incorporates in an analytically appealing way the effect of long-range 
dependence posed by the underlying nonlocal ordering. 

We start up with the conventional free energy expansion (Refs. [27] and [28]) 



F = F„+ / dx 



in vicinity of the critical point T ^ Tc- The term marked by Fn denotes 
the contribution from the normal phase. Equation (1) concerns with the real- 
space derivative V x"^ along the coordinate x. For the sake of simplicity, we 
shall also use the notation tp'^ for Vxip- The coefficient a = a{T — T^) in Eq. (1) 
is proportional to the deviation between T and Tc- Note that a — a{T — Tc) 
changes sign at criticality: This behavior does not depend on the nature of 
the symmetric and normal phases and mirrors the generic features of topology 
of the phase diagram. The system-specific information is contained in the 
parameters a and b. The value of a > in accordance with the fact that 
the symmetric phase occurs below Tc. The coefficient b > depends solely on 
the mass density of the material (but not on the thermodynamic temperature 
T). The bulk distribution of the order parameter extremizes the free energy 
functional in Eq. (1), leading to the classical Ginzburg-Landau equation [28] 

-—V\ilj + ailj + b\i/j\^ilj = 0. (2) 



The key issue about Eq. (2) is the infinitesimal — local — character of the 
ip^x) variation, manifest in the assumption that only a differential contribution 
oc |Va;'^|^ in the free energy density comes into play for T — > T^,. Statistically, 
this means that the blobs of the symmetric phase appear at random through- 
out the material as the temperature T approaches the critical range. The 
property of randomness is explicit from the convolution of Vxip with a Gaus- 
sian Gx{x) = exp[— a;^/A]/v^7rA, where A is the correlation (coarse-graining) 
length, physically corresponding to the typical size of the blobs: 

+ 00 

Vxi^*Gx{x) = J dytl^'^G^ix-y). (3) 

— oo 



At length scales x large compared to A, the Gaussian G\{x) in Eq. (3) can be 
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fairly approximated by the Dirac delta function 5{x), yielding Vx'ip*Gx{x) 
V xip * 5{x). The convolution Vxi^ * ^{x) is then exactly the local derivative 

Vx^ = ^'x- 

We now turn to the order parameter = "^{x) for the coexisting nonlocal 
symmetry. The interaction between ijj and $5 orders may be envisaged as a 
nonrandom appearance of the "superconducting" blobs for T ^ Tc over a 
broad range of scales a; ^ A. In this connection, the parameter Q=(a;) acquires 
the role of the Gaussian Gx{x) 5{x) in Eq. (3). Replacing S{x) by Q{x), 
one encounters the convolution Vx'4^*^{x), which substitutes the local deriva- 
tive Vx'ip * S{x) in the free energy expansion (1). We may now postulate the 
free energy expansion for a long-range correlated thermodynamical system at 
criticality in the generalized form 



F = Fn+ J dx A^\Vx^*'^\^ + a^\i)\^ + -hz 



(4) 



where $J = quantifies the underlying nonlocal symmetry. The coefficients 
a<3, and hc^ introduced in Eq. (4) replace /Am, a, and 6, respectively, in 
the conventional free energy expansion in Eq. (1). 

Our further interest is on the specific type of nonlocality consistent with a 
self-similar — fractal — organization. The focus on fractals [29] is motivated 
by a general tendency of complex systems to reveal, at or near a critical point, 
scale-invariant dynamical properties [16]. The phenomenon is often associated 
with the issue of self-organized criticality [30]. Self- organized critical behavior, 
due to multiscale Josephson coupling of the superconducting domains, was 
advocated for granular and polycrystalline superconductors in Ref. [31]. In 
the framework of our study, we assume the nonlocal symmetry has fractal 
support considered as a Cantor set [29] on the 1-dimensional Euclidean axis 
X. The fractal geometry of the support appears in the power-law behavior of 
the order parameter 9: This behavior, in turn, may be identified with the 
scahng of the two-point correlation function for the fractal distribution [29] : 

"^{x - y) = '^o\x - y\~^ (5) 



where Q'o is a normalization constant. The power exponent //. in Eq. (5) can 

further be expressed in terms of the Hausdorff fractal dimension of the 
Cantor set: 

11 = 1 -df. (6) 



By its definition [29] , the Hausdorff dimension of a Cantor set ranges from to 
1. In the latter case, the set occupies the entire Euchdean axis x. Accordingly, 



4 



the power exponent varies from a maximal value equal 1 to a minimal value 
equal 0. In view of Eq. (5), the convolution Vj;'0 * becomes 

+00 

V,V' * Q(x) ^% I dyi^'y \x - (7) 



The integration in Eq. (7) can be expressed in a suitable compact form by 
using the notion of a fractional derivative [32]. In fact, integrating by parts in 
Eq. (7), one gets 

+ 00 

* ^{x) ^% J dyiPiy) W^y\x - yl"''. (8) 

— oo 



The operation V-y applied to |a: — 1/| is equivalent with the derivative V^; over 
the parameter x: This derivative can then be taken out of the integral sign, 
yielding 

+ 00 

V,V^ * ^{x) = J dy^Piy) \x - y\-^'. (9) 

— oo 



Splitting the integration from —oo to +cx3 into two integrals, from —oo to x and 
from X to +00, and taking into account the reflection symmetry ip{—y) = ip{y) 
for the t/j order, from Eq. (9) one obtains 

X 

* ^(x) = 2^0 V, I dyi^iy) {x - y)'^. (10) 

— oo 



Setting the normalization 2^q — l/r(l — /i) in Eq. (5), we find * ^{x) = 
V^^t/j, where 

1 ^ 



is exactly the Riesz definition [32] of the fractional derivative of order < 
|U < 1, and r denotes the Euler gamma function. Note that V^^ is integro- 
differential operator for all < /i < 1. In the "integer" limit of — > 1, the 
operation in Eq. (11) is equivalent with the conventional first-order derivative 
Vx'- The proof rests on the Abel identities discussed in Ref. [33]. In terms of 
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fractional derivatives, the free energy expansion in Eq. (4) reads 

+00 



(12) 



Here we changed the subscript $5 to /i everywhere in Aq, cq, and Expression 
(12) leads to the issue of a fractional Ginzhurg- Landau equation, as we now 
proceed to show. In fact, varying the integral in Eq. (12) over the complex 
conjugate ■0* and considering and ■0* as independent order parameters, we 
have 



+00 



8F 



dx 



(13) 



Making use of the integration-by-part formula 

+00 +00 

j dy cpi{y)Vyip2{y) = j dy ip2{y)V%(pi{y) 



(14) 



from Eq. (13) one arrives at 



-|-oo 



(15) 



yielding, in view of the extremum 6F — 0, 



(16) 



Varying the integral in Eq. (12) over t/j leads to the conjugate equation 

A,v^_,w,r + a^r + = o (ir) 



which is physically identical to Eq. (16). Equation (16) can be considered as 
a fractional generahzation of the Ginzburg-Landau Eq. (2). The fractional 
Ginzburg-Landau Eq. (16) determines the bulk distribution of the order pa- 
rameter ip in the presence of a coexisting nonlocal symmetry, whose support 
is a Cantor set of the HausdorfT dimension dj = 1 — fi. The coefficient 
consistent with the fractional diffeo-integration in Eqs. (12) and (16) could be 
defined by 



4m 



(18) 
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where A is the microscopic correlation length. In the conventional case of 
// — > 1, the value of reduces to h^/Am. 

Let us now discuss a likely form for the coefficient a^u in Eq. (16). Without loss 
of generality, we may admit that is a linear function on thermodynamic 

temperature T. (We assume that all explicit nonlinearities have already been 
allocated to a single term traced by 6^.) Let cross zero at some T — 
i.e., 



Equation (19) is a simplest extension of the conventional formula a — a{T—Tc) 
central to the classical Landau theory (Refs. [27] and [28]). An important is- 
sue about Eq. (19) is that the cross-over temperature where changes 
sign may not coincide with the transition temperature when the effects of 
nonlocality play a role (i.e., when the index /x is smaller than 1). From the 
standpoint of a formal derivation, the value of Tc enabling the occurrence of 
a nontrivial ip order in the presence of the competing nonlocal symmetry de- 
couples from owing to the fractional nature of the Riesz operators V'^_^ 
and V^^ which integrate the response from the heavy-tailed correlation func- 
tion '^{x — y) in Eqs. (11) and (16) and thereby contribute (along with the 
parameter a^) into an amplitude in front of ■0. (This amplitude then balances 
the nonlinear term fe/^lV^P below the transition point.) Physically, the devia- 
tion between and has the sense of a characteristic energy of correlations 
contained in the fractional Laplacian, ^4^ V'^_^V'^^. By order of magnitude. 



where F is the Euler gamma function. As ji tends to 1, the deviation in Eq. 
(20) vanishes: F(l — /x) — > oo leading to — > T^. The conventional expression 
a — a{T — Tc) is then recovered from Eq. (19). In this limit, the nonlocal 
symmetry confines on a set whose Hausdorff measure is zero, i.e., df for 
/X — > 1. This behavior reproduces the assumptions of the classical theory [28]. 

In the opposite limit of // — > 0, the fractal support underlying the nonlocal 
symmetry extends to the whole of the 1-dimcnsional coordinate space whose 
Hausdorff measure is now maximized to he df 1. This almost regular distri- 
bution bears features enabling to associate it with a "stripe" [25], a long-range 
ordering that seems to mediate superconductivity in complex materials such 
as copper-oxide compounds and their derivatives [26]. 

As /i — s> 0, the gamma function in Eq. (20) drops to a minimal value equal 
1, enabling a considerable deviation between Tc and T^. It can be conjec- 
tured that, in general, this deviation accounts for the "anomalous" transition 



an = Oi{T -Tfj^). 



(19) 



Tc-T,r.A,/aX''^r\l-f,) 



(20) 
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temperatures observed in high-temperature superconductors. A self-consistent 
estimate of the parameter may then be obtained from Eq. (20). A key step 
is the expression (Refs. [28] and [34]) 



which relates the quantity a to the system-specific characteristics such as 
electron energy at the Fermi edge, ep (typically, of the order of few eV). The 
coefficient T in Eq. (21) depends on the nature of the microscopic pairing 
mechanism (e.g., phonon, fracton, or exciton) (see Refs. [18] and [35]). In the 
case of phonons — acoustic or (quasi) acoustic modes serving as an interface 
for the conventional Bardeen-Cooper-Schrieffer (BCS) picture of superconduc- 
tivity in regular crystals [36] — the value of T ~ 7 (Refs. [28] and [34]), but 
it can be much smaller (of the order of 1 or even less) in the Little's exciton 
scenario, in which phonons are replaced by other electrons [37]. Combining 
Eqs. (20) and (21), we find, in the extreme of — > 0, 



provided that ^ T^. Assuming, further, the correlation length A is of 
nanoscales — the range typically found for the high-temperature superconduc- 
tors [38] — from Eq. (22) one concludes a characteristic transition temperature 
Tc ~ 100 K where the estimate T ~ 7 has been used. Non-BCS pairing pro- 
cesses corresponding to smaller values of T would imply a considerably higher 
critical temperature T^, up to room-temperature values. A feasible role of such 
processes in two-dimensional electron systems with multi-scale long-range cor- 
relations is addressed in Refs. [35] and [39]. A fractional kinetics constituent 
in the occurrence of the ensuing superconducting fluid state will be the sub- 
ject of a forthcoming publication. We believe this constituent is crucial for 
the unconventional superconductivity in materials with complex microscopic 
organization. 
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ferent aspects of fractional kinetics. This study was sponsored by the Science 
Support Foundation, INTAS Grant 03-51-3738, RFBR project 03-02-16967, 
and "Scientific School" Grant 1739.2003.2. 
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